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Abstract. The chiral symmetry breaking properties of the Tayler instability are discussed. 
Effective amplitude equations are determined in one case. This model has three free parameters 
that are determined numerically. Comparison with chiral symmetry breaking in biochemistry is 
made. 
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1. Introduction 

An important ingredient to the solar dynamo is the a effect. Mathematically speaking 
a is a pseudo scalar that can be constructed using gravity g (a polar vector) and angular 
velocity H (an axial vector): g • is thus a pseudo scalar and is proportional to cos0, 
where Q is the colatitude. This pseudo scalar changes sign at the equator. This explanation 
for large-scale astrophysical dynamos works well and therefore one used to think that 
the existence of the a effect in dynamo theory requires always the existence of a pseudo 
scalar in the problem. This has indeed been general wisdom, although it has rarely been 
emphasized in the literature. That this is actually not the case has only recently been 
emphasized and demonstrated. One example is the magnetic buoyancy instability in 
the absence of rotation, but with a horizontal magnetic field B and vertical gravity fj 



being perpendicular to each other, so the pseudo scalar g ■ B vanishes (jChatteriee et al 



201 If). Another exa mple is the Tayler instability of a purely toroidal field in a cylinder 
(jGellert et al. 2011 1. Thus, the magnetic field is again perpendicular to all possible polar 
vectors that can be constructed, for example the gradient of the magnetic energy density 
which points in the radial direction. In both cases, kinetic helicity and a finite a, both 
of either sign, emerge in the nonlinear stage of the instability. In the former case, the 
a tensor has been computed using the test-field method. In the latter, t he components 



of th e a tensor have been computed using the imposed- field method (see lHubbard et al 
20091 for a discussion of possible pitfalls in the nonlinear case). 



The purpose of the present paper is to examine spontaneous chiral symmetry break- 
ing in the Tayler instability and to estimate numerically the coefficients governing the 
underlying amplitude equations. This allows us then to make contact with a model sys- 
tem of chemical reactions that can give rise to the same type of spontaneous symmetry 
breaking. 

The connection with chemical systems is of interest b ecause t h e que stion of sponta- 
neous symmetry breaking has a long history ever since iPasteun (|l853l ) discovered the 
preferential handedness of certain organic molecules. The preferential handedness of 
biomolecules is bel ieved to be the result of a bifurcation event that took place at th e 
origin of hfe itself ( Kondepudi fc Nelson 1984 : Sandars 2003 : Brandenburg et al. 2005f) . 
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Figure 1. Evolution of magnetic helicity for two initial conditions, differing only in the parity 
of their initial perturbations. After the exponential growth magnetic helicity levels off. In the 
inset a detail of the exponential growth phase. Here, R = Sin is used. 

2. Numerical simulations 

Our setup consists of an isothermal cylinder with a radial extent from Sin to Sout and 
vertical size h. We solve the time dependent ideal MHD equations with periodic boundary 
conditions in z, reflection in s and periodic in if and a resolution ranging from 64'^ to 
128^ in the three directions. 

The azimuthal field in the basic state is taken of the form 

= bo (s/so)exp[-(s - so)^/o-^] 

with bo being a normalization constant; the axial field Bz is chosen to be zero. In the basic 
state, the Lorentz force is balanced with a gradient of pressure, and we have checked that 
our setup was numerically stable if no perturbation was introduced in the system. For 
the actual calculations we have chosen h = 2, Sin = 1, Sout = 3, sq = 2 and = 0.2. The 
sound speed is assumed to be much larger than the Alfen speed (w ten times), similar 
to what happens in stellar interiors. 

At the beginning of the simulation we perturb the magnetic field. We add a pertur- 
bation of amplitude 10~^ that of the background field. The perturbing field has a given 
helicity that is either positive or negative. During the development of the instability we 
observe a net increase of the helicity, as shown in Fig. [1] where we plot time series of 
the normalized magnetic helicity, which exhibits an initial exponential growth, reaches a 
peak and then levels off. 




3. Amplitude equations 

The linear stability analysis of this instability shows that there exists helical growing 
modes. But the left handed and right handed modes have exactly the same growth 
rate independent of their helicity. Hence the growth of helical perturbations cannot be 
described by a linear theory. However a weakly nonlinear theory is able to describe it as 
we show below. Let us begin by considering two helical modes of right handed and left 
handed variety respectively each of which satisfy the Beltrami relation V x i? = AR and 
V X L = —AL . We can deal with the Fourier transform of these modes, given by 

L{x) = J L{q)d'^q and R{x) = j R{q)d'^q (3.1) 

For the left helical mode, total helicity and energy are given by 



El 

2 



^ J L^{x)d'^x ^^Jl- L*d'^q and TiL = J h ■ V x Ld'^x =^ -2AEl, (3.2) 
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where * denotes complex conjugation. We then have E = +Efj being the total energy 
and H = Hl + Hr the total helicity. An analogous relation holds for the right-handed 
helical mode too. 

In the weakly nonlinear regime the evolution of these modes can be described by 
general equations of the form: 



'dt 



(3.3) 



5C_ dR _ 5C 

Jl 'dt ^Jr' 

where the Lagrangian C can often by written down from symmetry considerations. In 
the present case one has to consider the fact that under parity transformation L can 
R interchanges into each other. With this a dditional symmetry the simplest Lagrangian 
takes the following form (jFauve et al.lll99ir) 



C[L, R] = 



7 



\R\- 



— fi \L 



\R[ 



fi,\L\^\R\^d'^q\ d'^q, 



(3.4) 



The coefficients 7, /i and /x* cannot be found from symmetry considerations. Note that 
in order to show the simplest form, in writing down the Lagrangian we have ignored 
dissipation. This gives rise to the following set of amplitude equations, 



8T 



\ dR / 



R- 



R. 



(3.5) 



For certain range of paramet ers these coupled equations allow the growth of one mode 
at the expense of the other (IFauve et al.l 1199111 . a phenomenon known to biologists by 
the name "mutual antagonism" (|Franklll953f l 

Using Eqs. p.2|) and p.5|) and defining H = H/2A we can obtain evolution equations 
for E and H as 



— = 2-/H - AfiEH. 
dt 



(3.6) 
(3.7) 



Hence, by calculating the total energy and helicity from direct numerical simulations 
(DNS) we can determine the unknown coefficients 7,^ and /z*. 

To determine the coefficients 7, fi, and /i*, we define the instantaneous logarithmic 
time derivatives of E and H, 7^ = \d\aE/dt and 7/f = idlniJ/df, so we have 

7 = 7^ + 2iiE, M = (7 - 1h)/2E, fi, = [(7 - 1e)E - fi{E^ + H^)]/{E^ - H^). (3.8) 

The result is shown in Fig. [2j where we can identify first the value of 7 « 14 during the 
initial linear growth phase of the instability, and then the values /i « 10 and fi, ^ 7 
during the nonlinear stage. 



4. Conclusions 

The present work has demonstrated that the Tayler instability can produce parity- 
breaking and that it is possible to determine empirical fit parameters that reproduce 
the nonlinear evolution of energy and helicity. So far, no rigorous derivation of the am- 
plitude equations exists, so this would be an important next step. Comparing with the 
chiral symmetry breaking instabilities in biochemistry, an important difference is that in 
the present equations the nonlinearity is always cubic, while in biochemistry the domi- 
nant nonlinearity tends to be quadratic. In this light, it would be useful to assess more 
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Figure 2. Time dependence of 7, /x, and /is, normalized in terms inner radius and sound 
speed. The red lines give the fit results 7 « 14, /i ~ 10, and /x* « 7 in the appropriate units. 



closely the possible differences between biochemical and magnetohydrodynamical sym- 
metry breaking. 
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